Leading Temperature Corrections to Fermi Liquid Theory in Two Dimensions 
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We calculate the basic parameters of the Fermi Liquid: the scattering vertex, the Landau in- 
teraction function, the effective mass, and physical susceptibilities for a model of two-dimensional 
(2D) fermions with a short ranged interaction at non-zero temperature. The leading temperature 
dependences of the spin components of the scattering vertex, the Landau function, and the spin 
susceptibility are found to be linear. T-linear terms in the effective mass and in the "charge-sector" - 
quantities are found to cancel to second order in the interaction, but the cancellation is argued not 
to be generic. The connection with previous studies of the 2D Fermi-Liquid parameters is discussed. 
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The question of the low-energy behavior of two- 
dimensional (2D) Fermi Liquid (FL) is of long-standing 
and of fundamental importance. One important motiva- 
tion has been the non-Fermi-Liquid behavior observed in 
high-T c superconductors above T C .EI Interest has grown in 
recent years also because these leading corrections pro- 
vide the "bare" temperature dependence of the param- 
eters in theories describing quantum critical phenomena 
in metals!] Indeed a numbgcof surprising experimental 
resultstrQ have been arguedtm to imply an unusal under- 
lying temperature, momentum or frequency dependence 
of electronic susceptibilities. 

Rather surprisingly, the issue of the leading tempera- 
ture corrections to Fermi Liquid Theory (FLT) remains 
controversial at certain points. For example, the text- 
book Sommcrfcld expansion suggests that physical quan- 
tities may in general be expanded in powers of T 2 . How- 
ever, this is known not to be true. In particular, it was 
foundEI that the leading temperature correction to the 
specific heat coefficient 7 = C/T was T 2 InT in d=3 spa- 
tial dimensions and T in d=2j9. Whether the spin and 
charge susceptibilities display similarly anomalous (i.e., 
non-T 2 ) temperature dependence is a subject of a con- 
tradictory literature: see, e.g., Ref. [ p|, discussion and 
references there in. For a most recent reassesment of such 
results see Ref. [ . The prevailing conclusion was that 
of Carneiro and Pethicklj who found no leading T 2 In T 
correction to the spin susceptibility of the 3D FL. Their 
arguments imply that terms oc T are absent in 2D. 

The heuristic argument for the absence of anomalous 
terms in T or in q in response functions is that although 
these terms are known to occur in individual diagrams, 
they cancel in physical quantities due to Ward identities. 
We note that in the existing literature on this point it 
is assumed that the crucial coupling is between quasi- 
particles and long-wavelength collective modes. How- 
ever the possibility of ll 2kp singularities", i.e., anomalous 
temperature terms coming from processes involving large 
(~ 2fcp) momentum transfers, have been discussed in the 
context of semiconductor physics. Stern was the first to 
noteO that in a 2D electron gas the electron scattering 
rate was proportional to T due to 2kp effects. The conse- 



quences of the 2kp effects for the leading T-dependence 
of 2D FL quantities seem not to have been considered in 
the literature. In this paper we present an analysis taking 
into account both 2kp effects and Ward identities. 

The issue of the leading correction to FLT have re- 
cently been revived-by two different papers. Belitz, Kirk- 
patrick and VojtaEJ presented perturbative calculations, 
mode-coupling arguments and power counting estimates 
which showed that the leading q dependence of the spin 
susceptibility (but not the charge susceptibility) was \q\ 
in 2D and q 2 Inq in 3D. They did not find the analogous 
T-correction explicitly, but concluded however that one 
should generally expect a linear T-term in the 2D FL 
susceptibility (T 2 lnT in 3D). This dependence has im- 
portant implications for the theory of the quantum crit- 
ical metallic ferromagnet.E3 They also focused on long- 
wavelength contributions. _ 

In the other study, Senechal and one of usE3 predicted 
the occurence of the linear T-corrections to the FL ver- 
tices from one-loop Renormalization Group (RG) calcu- 
lations based on a 2D effective action. However, the be- 
havior of other FL quantities was not determined. 

To elucidate the question in the most transparent way, 
we apply perturbation theory for 2D contact-interacting 
spin-i fermions, starting from a microscopic action. Al- 
though Landau FLT is not a perturbative theory, for suf- 
ficiently weak repulsive interactions one should be able 
to find the parameters of a stablepEL in terms of cou- 
pling series. Since both papersEjuJ predict the linear 
T-terms appearing in the second order of the effective in- 
teraction, this effect should be seen perturbatively at sec- 
ond order in the miscroscopic interaction coupling. We 
present what is apparently the first exact calculation of 
the leading T-dependence of the effective mass, Landau 
parameters and response functions of a 2D electron gas, 
to second order in the interaction strength, including all 
channels and all momentum processes. 

• The model: We treat interacting fermions at finite 
temperature, in the standard path integral Grassmannian 
formalism.Eil The partition function is given by the path 
integral Z = J VipV^p exp(So + S- in t), where 
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S = Va(l)[i«i+/i-e(ki)]Va(l) (1) 
J(l) 

We have adopted the condensed notations: (i) = (kj,cjj) 
and = jj J Y^Ui > P 1S the inverse temperature, fx 
the chemical potential, loi the fermion Matsubara fre- 
quencies and tp a (X) a two-component Grassmann field 
with a spin index a. Summation over repeated indices 
is implicit throughout this paper. We set fcg = 1 and 
% = 1. We consider mainly electrons with the bare spec- 
trum of a free gas e(k) = k 2 /2m and the circular Fermi 
surface, but discuss the consequences of generic spectra 
and a non-circular Fermi surface below. We take 
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V; a (l)#(2)V> 7 (3)V> e (4)T 7 ^ 



<S (2+1) (l + 2-3-4) 



(2) 



Here u > corresponds to repulsion, vq = m/2ir is the 
free 2D density of states per spin, the spin antisymmetric 
operator T"f = 5 ae 5p 



— 5 ai 8j3 e . For this model we 
calculate the vertices, effective mass, and the uniform 
response functions. 

• Four-point 1PI vertex and the FL vertices: 

The 1PI vertex T(l, 2; Q) is defined in the standard way 
(see, e.g., Ref.[ Q), where the transfer vector Q = (q, fi), 
and Vt is a bosonic Matsubara frequency. To shorten 
notations we denote operators in spin space with a cir- 
cumflex. To find the FL vertices we need to calculate 
r(l, 2; Q) in the limit of zero transfer Q with incoming 
momenta lying on the Fermi surface. Since the Fermi 
surface is circular, the vertex can be parametrized by 
the relative angle between incoming momenta. It is 
knowr£3'ta that the limit Q — -> is not unique. We de- 
fine two vertices f«(0i 2 ) = r(0i 2 ;q -> 0, fl = 0) and 
f °(0 12 ) ee f (0 12 ; q = 0, il -» 0), which can be related to 
the components of the physical scattering vertex (A, B) 
and the Landau interaction function (F, G), respectively. 
Namely:E3 



(3) 



where Vr = m* /2w, Z is the field renormalization con- 
stant, a a are Pauli's matrices. Two components (F, G) 
of the Landau function are defined by analogous equa- 
tion, with the substitution q i— » Q, A i— > F, B i— » G. 

• Scattering vertex and Landau function: The 
one-loop approximation for r(l,2;Q) in diagrammatic 
form is given in Fig. 1. In this approximation we calcu- 
late the FL vertices (scattering vertex, Landau function) 
using definition d). At the one-loop level we can put 
V R = v q ( m * = m ) and Z = 1. By doing the direct 
analytic evaluation of the each diagram's contribution to 
the vertex we find the Fourier components of the angular- 
dependent FL vertices in terms of the temperature series. 
This series comes from the contributions of the ZS' and 
BCS Ipops. The details will be given in a companion 
paperEO We find the leading temperature corrections to 
first two Fourier components of the vertices: 



6A Q = SF = 8 A x = SFx 



5B = SG = -SBx = -SGx 
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(4) 
(5) 



We should mention that taken separately, each contribu- 
tion of the ZS'- or BCS bubble gives a leading T-term 
to the Fourier components of the vertices. However, a 
cancellation of such terms coming from two graphs oc- 
curs in the "charge sector" (i.e., in A, F components), 
while the linear T-terms survive in the "spin sector" (B, 
G components). The temperature dependence of the ZS' 
contribution to the Fourier components of the vertices 
comes from integration around the "effective transfer" 
through the loop |ki — k 2 | ~ 2kp, i.e., when incoming 
momenta ki ~ — k 2 . In the same vein, the tempera- 
ture dependence of the BCS contribution into the Fourier 
components of the vertices comes from regions of small 
ki + k 2 , i.e., again when ki ~ — k 2 . In other words, 
the temperature dependence comes from what we called 
previously u 2kp-eSects" . 

We expect that the cancellation of the T-linear terms 
in the charge sector of the vertices (f|) is an artefact of 
our simple model calculation in which all three one-loop 
terms have the same factor u in front of the bubble con- 
tributions. Had we had a bare coupling function of, say, 
two incoming momenta and transfer, then the coupling 
factors would have been different in each of the three 
graphs, and, the anomalous T-corrections would not have 
cancelled. 

The linearity of the leading T-corrections to the ver- 
tices seems to be generic. The same T-dependence 
(apart from presumably model-sensitii££ prefactors) was 
obtained in the previous RG analysisEj of the effective 
action for 2D spinlcss fcrmions with a linearized one- 
particle spectrum and a momentum-dependent coupling 
function. According to Ref. [ [T(| such temperature be- 
havior can be understood from dimensional arguments. 

• Effective mass: It is defined by the following equa- 
tion 
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To second order we have 
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+ eV) (7) 
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Using then two Ward identities followingJijom the charge 
conservation and Galilean invariance,t3ll3 the above 
equation can be written as: 
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r^(l,2;fi^0)A(2) 
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where 
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4cosh 2 (i^ k J 



(9) 



Within our accuracy we can use the one-loop approxi- 
mation for the vertex in Eq. (^). One can easily verify 
that in the the zero-temperatute-, limit Eq. (j^) recovers 
the standard result of the FLT,EJ i.e., m*(T = 0)/m = 
1 + -Fx (T = 0). A straightforward extension of this rela- 
tionship to finite temperatures like m*(T)/m = l+Fi(T) 
is not valid since according to Eq. (|J) m* (T) contains an 
extra contribution from the "off-shell" integration over 
&2 (£k 2 ) normal to the Fermi surface, albeit the factor 
(3/ cosh 2 (/3£k 2 /2) makes this contribution well-localized 
near the Fermi surface. In other words the vertex en- 
tering the r.h.s. of Eq. (||) is not exactly the FL ver- 
tex F(T) (up to the normalization factor) as it is de- 
fined in the FLT, since one of its momenta (namely, 
is not confined to the Fermi surface. After calculations 
we find that the linear-temperature terms, coming essen- 
tially from two one-loop contributions (ZS', BCS) to the 
vertex, cancel, resulting in 



m* , In 

m 2 



0{u 2 T 2 



(10) 



In a close analogy with the cancellation of the linear- 
temperature terms in the Fourier components of the FL 
vertices A (T), here the cancellation occurs between ad- 
ditive linear-T corrections coming from both "on-shell" 
(i.e., linear T-term coming from the 2/c^-contribution 
to the vertex) and "off-shell" (i.e., the small- momentum 
contribution) integrations in two diagrams. Moreover, 
the "on-shell" ("off-shell") T-term of the ZS' graph can- 
cells the "on-shell" ("off-shell") T-term of the BCS graph, 
correspondingly. We expect the cancellation does not oc- 
cur at higher orders in the interaction. We have also 
evaluated (]?]) for a generic 2D Fermi surface without ex=. 
plicitly using Ward identities, finding a T-linear term.EZl 
The result may be expressed as the sum of two terms, one 
arising from 2k p processes and one from long- wavelength 
processes. The two contributions cancel for a circular 
Firmi surface, but not generically. 



x FLT = ^ jj^- , expanded up to the third order over the 
interaction. Add ing into consideration the Ward identity 
following from the total spin conservation, we derived for 
the uniform spin susceptibility (for details see Ref. [ |l7|): 



X = — — <1 + u 



4tt 



[l + u 2 + h -g } 



(12) 



5o = ~ [ A(l)a>^(l, 2; Q - 0)A(2) 

where g stands for the gyromagnetic ratio. Once again, 
one can see that in the zero-temperature limit the above 



equation gives X{T = 0) = ^- (1 + G 2 , + F 1 - G ) repro- 
ducinj 

resultl 



thus the second-order expansion of the the FLT 

vFLT _ mg 2 1+F\ 
^ 4tt 1 + Go ' 

We were able to analytically calculate the integrals on 
the r.h.s. of Eqs. jll|Jl^ ) in the leading order of their 
temperature dependence. We found that the leading lin- 
ear T-corrections, which can be traced back to the ZS'- 
and BCS-loop contributions to the vertex, cancel in each 
of the integral terms /o and /i in Eq. (|lTJ) separately. 
The result for the density response is: 

k=-{1-u-\u 2 + u 2 \n^) + 0(u 2 T 2 ) (13) 

7T 2 kp 

where A 3> kp is the ultraviolet cutoff we intpaduced to 
regularize the BCS loop. Let us finally remindllS that the 
compressibility K — x/n 2 , where n is electron density. 

We may calculate the spin susceptibility in the same 
way. In this case the second integral term g in Eq. (|l^ ) 
does not contain the contribution of the ZS' loop, so the 
linear T-term coming from the BCS loop survives. Thus 
the spin susceptibility has a linear leading temperature 
correction: 



X{T) 



T 
~Ef 



(14) 



where X(0) = 2S-[1 + u - u 2 {\n |A - §)], and the first 



Calculations of the order u 2 term in the free energy omitted term is 0{u T ) 



show that analogous cancellations occur and there is no 
T-linear term in the specific heat coefficient 7, contrary 
to the results of Ref. [§. 

• Response functions: Using the same Ward iden- 
tities as in the effective mass calculation, we found for 
the dynamic zero-transfer limit (fi = 0, q — > 0) of the 
density response function: 



It is useful to kee 
functions in Eqs. (|l 



x= ^{l + u 2 + /l _ /o j 



(11) 



7r f kik 2 
in ./(i,2) k l 



/, / ^A(l)r2g(l,2;fi->0)A(2) 



fo = ~ [ A(l)r$(l,2;n-0)A(2) 

m J(l,2) 

At T = we can read off from Eq. @ x(T = 0) ^ 
^(1 + T 2 +Ti-T ) which is nothing but the FLT result^ 



in mind that albeit the response 
are explicitly expressed in terms 
of the vertex only, those contributions indeed entangle 
both "purely vertex" corrections and self-energy correc- 
tions. The latter are just expressed in terms of the vertex 
via the Ward identities. 

Let us come back to the argument for the cancellation 
of temperature terms in the response functions, appeal- 
ing to the Ward identities. These identities should work 
at each level in terms of the coupling expansion. We 
have calculated the vertices at the one-loop level 0(u 2 ). 
Through the Ward identities the self-energy correction 
were taken into account with the same accuracy. There 
are no more terms of the order 0(u 2 ) to cancel the tem- 
perature dependence (14). Thus, the linear temperature 
dependence of susceptibility (or weaker temperature de- 
pendence of the compressibility) does not contradict the 
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exact Ward identities known to us, moreover in our re- 
sults for the response functions both vertex and self- 
energy corrections are included on the same footing by 
using the Ward identities. 

• Conclusions: In this paper we have systematically 
examined the leading temperature corrections to FLT in 
two spatial dimensions. Our results reveal the crucial 
importance of 2kp processes mainly neglected by other 
workers. We find for a 2D Galilean-invariant FL that to 
order u 2 the leading T-dependence of the parameters in 
the spin sector is T, for the others it is T 2 . 

We find that the standard relationship between the 
effective mass and the Landau parameter in a Galilean- 
invariant FL is violated by finite-temperature terms. The 
structure of the formulas shows that those terms arise 
from states very near to the Fermi level ( "off-shell" quan- 
tities) beyond the scope of the classical FLT derivations. 

The particularly interesting new result we found is the 
leading linear temperature dependence of the spin sus- 
ceptibility ([l4]). According to the perturbative calcula- 
tions of Bclitz et a^Ej the 2D FL susceptibility has a 
leading linear correction in |q| at T = with a positive 
coefficient which is of the second order in interaction, 
i.e., their result has a structure of Eq. (|l4|). Comparison 
between our Eq. ( |i~4| ) and the prediction of Ref . [ |l0| for 
<T(q, T = 0) meets heuristic expectations of a reciprocity 
between small T- and q-dependencies. 

For more realistic models of electrons in (quasi)-2D 
crystals, i.e., for various tight-binding spectra and fill- 
ings, the free-gas-like square-root 2kp singularities (with 
kp depending on a chosen direction in_q-space) are known 
to exist in the Lindhard functions.E£l We think this is 
enough to result in linear T-terms in physical quantities 
analogous to what we found in this study. We argue 
that the rather accidential cancellation of the T-terms 
in some FL parameters is special to second order per- 
turbation theory and a circular Fermi surface, while the 
leading linear temperature corrections are a generic fea- 
ture of the 2D FL. 

We hope our results may be experimentally tested in 
real 2D FL systems. For example, a very naive fit of 
the temperaturejdependence of the spin susceptibility in 
Si'2Ru04 systemtLj when it is in the 2D metallic regime 
(above 3D crossover temperature) shows that the data 
are compatible with the form (jl4|). We hope our re- 
sults stimulate a more detailed examination of the lead- 
ing temperature dependences of response functions in 2D 
systems. 
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FIG. 1. Diagrammatic equation for the four-point vertex 
at one-loop level. The one-loop graphs are called ZS, ZS', and 
BCS in the order they appear on the r.h.s. of this equation. 
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